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ABSTRACT 

High-sensitivity measurements of the linearly-polarized solar limb spectrum produced by scattering pro- 
cesses in quiet regions of the solar atmosphere showed that the Q/I profile of the lithium doublet at 6708 A 
has an amplitude ^10""* and a curious three-peak structure, qualitatively similar to that found and confirmed by 
many observers in the Na i D2 line. Given that a precise measurement of the scattering polarization profile of 
the lithium doublet lies at the limit of the present observational possibilities, it is worthwhile to clarify the phys- 
ical origin of the observed polarization, its diagnostic potential and what kind of Q/I shapes can be expected 
from theory. To this end, we have applied the quantum theory of atomic level polarization taking into account 
the hyperfine structure of the two stable isotopes of lithium, as well as the Hanle effect of a microturbulent 
magnetic field of arbitrary strength. We find that quantum interferences between the sublevels pertaining to the 
upper levels of the D2 and Di line transitions of lithium do not cause any observable effect on the emergent 
Q/I profile. Our theoretical calculations show that only two Q/I peaks can be expected, with the strongest 
one caused by the D2 line of ^Li i and the weakest one due to the D2 line of ''Li i. Interestingly, we find that 
these two peaks in the theoretical Q/I profile stand out clearly only when the kinetic temperature of the thin 
atmospheric region that produces the emergent spectral line radiation is lower than 4000 K. The fact that such 
region is located around a height of 200 km in standard semi-empirical models, where the kinetic temperature 
is about 5000 K, leads us to suggest that the most likely Q/I profile produced by the sun in the lithium doublet 
should be slightly asymmetric and dominated by the ^Li i peak. 
Subject headings: Polarization - Scattering - Sun: atmosphere 
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1. INTRODUCTION 

One of the interestin g peculiarities of t h e linea rly-polarized 
spectrum observed by IStenflo & Relief (11997b in quiet re- 
gions close to the solar limb is that a variety of spectral 
lines from minority species, which produce almost negligi- 
ble absorption features in the Fraunhofer spectrum, nonethe- 
less stand up with significant contrast in fractional linear po- 
larization (i.e., in the Q(A)/I{A) spectrum)'. For example, 
molecules make a significant contribution to the structural 
richness of this so-called Second Solar Spectrum, and it is 
of interest to note that their magnetic sensitivity through the 
Hanle effect has facilitated the exploratio n of the sun's hidden 
magnetism (see the detailed review by ITruiillo Bueno et al] 
l200a and more references therein). Spectral lines from rare 
earths, such as those from ionized cerium, are also of inter- 
est for obtaining information on unresolved, tangled magnetic 
fields in the "quiet" solar atmosphere (see Manso Sainz et al. 
I2OO6I) . The main aim of this paper is to investigate the phys- 
ical origin and the diagnostic potential of what is proba- 
bly the weakest Q/I spectral feature produced by a minor- 
ity species in t he whol e Second Solar Spectrum: that ob- 
served by Sten flo et al.l (j2000) in the Li i resonance dou- 
blet at 6708 A (see Figure [1]). Although it is generally ac- 
cepted that the Second Solar Spectrum is due to radiatively 
induced population imbalances and quantum coherences in 
the atoms and molecules of the solar atmosphere, and sev- 

' Note that the first Stokes parameter, is the specific intensity at a 
given wavelength, while Stokes Q(A) represents here the intensity difference 
between linear polarization parallel and perpendicular to the closest solar 
limb. 



eral interesting spectral lines have been successfully modeled 
(e.g., Manso Sainz & Trujillo Bueno 2003), much remains to 
be done to fully understand the Q/I profiles observed in many 
other spectral lines. 

As seen in Figure [T] th e Une center value of the observed 
Q/I profile published by Stenflo et alJ (120001) is only 2x10 
and its detection required to sacrifice completely the spatio- 
temporal resolution in order to be able to push down the nns 
noise to below the 10"^ level. The FWHM of the observed 
Q/I profile is about 190 mA. The most peculiar feature of 
the observed signal is its three-peak structure, which bears 
a qualitative resemblance to the triple peak structure of the 
Q/I profile observed in the Na i D2 line. Like the sodium 
doublet, the Li i one at 6708 A results from a D2 type tran- 
sition (Ju - 3/2 and J( = 1/2) and a Di type transition 
(Ju = 1/2 and J( = 1/2), with the difference that lithium 
has two isotopes: ^Li (with nuclear spin / = 3/2 and a rel- 
ative meteoritic abundance of 92.41%) and ''Li (with 7=1 
and a relative meteoritic abundance of 7.59%). The fact that 
quantum interferences between the J„ -3/2 and 7„ = 1/2 
levels are indeed important for understanding the Q/I pat- 
tern ob served at wavelength s around the sodium D2 and Di 
Unes (Stenflo & Kelleil ll997l;lLandi DeglTnnocentil 19981) led 
Stenflo et al. (2000) to argue that the Q/I feature shown in 
Figure[T]must be seriously influenced by a quantum mechani- 
cal superposition of the scattering transitions in the D2 and Di 
lines of lithium. 

An accurate measurement of the very weak linear polar- 
ization signal of the Li i doublet at 6708 A is so difficult 
that we believe that it is very important to investigate care- 
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TABLE 1 

Physical properties of the two stable isotopes of lithium 



Isotope Abund. / Isotope Shifts (MHz)° HPS Constants (MHz) 

(%f Di D2 ^51/2 ^^3/2 

s_ 

*Li 7.59 1 -10533. IS"-' -10534.93''' 152.1368393'' 17.375= -1.155= 

-0.010= 

'Li 92.41 3/2 reference isotope 401.7520433'' 45.914' -3.055^ 

-0.22 if 



°A positive i s otope s hift means that the line is shifted to higher f requencies with r espect to the reference isotope. iRalchenko et alj ( 120081) ; 
iScherf et al.l ( fl99^ ; - iBeckmami et al.l ( fT974l) ; lOrth etal] ( fl974l) ; lOrth et al.l ( fl97^ . 
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Fig. 1 . — Intensity spectrum (upper panel) and fractional linear polarization 
(lower panel) of the Li i 6708 A doublet, as observed by Stenflo et al. (2000). 
The wavelength positions of the two components of this doublet (which are 
of the D2 and Di type), due to the most abundant lithium isotope ('Li), are 
indicated by the dashed and dotted lines, respectively. The r ecording was 
made on 1996 September 8. near the south polar limb. From [Stenflo et all 
<2000l) . 

fully w hether or not the 2/^ profile observed bv lStenflo et al.l 
(|2000) can be confirmed by theory. We think that this can 
be achieved with a particularly high level of confidence for 
this lithium doublet, because its extreme weakness in the in- 
tensity spectrum suggests that radiative transfer effects play 
no significant role on the shape of the Q/I profile. In- 
deed, the anisotropic radiation pumping is basically that due 
to the solar continuum radiation. Moreover, since the in- 
cident radiation field is practically flat across the lithium 
doublet, we can rigorously apply the complete redistribution 
theory of spectral Une polarization (se e the monograph by 
iLandi Degl'Innocenti & Landolfil 12004'. hereafter referred to 
as LL04), either neglecting or accounting for quantum inter- 
ferences between the hyperfine structure (HFS) F-levels per- 
taining to the Ju = 3/2 and 7„ = 1/2 levels. 

As we shall see, the application of the above-mentioned 
quantum theory of spectral line polarization leads to the con- 
clusion that only two peaks can be expected for the Q/I pro- 
file of the Li I resonance doublet at 6708 A. Interestingly, we 
find that the strongest peak at 6707.75 A is caused by the D2 
line of ^Li i, while the weakest one at about 6707.9 A is due 



to the D2 line of *Li i. Since in the sun the line opacity of 
this lithium doublet is negligible with respect to that of the 
continuum, contrary to the case of strong lines like Na i Di 
and D2, or Ca 11 H and K, these lithium lines neither show 
extended wings in the intensity spectrum nor in the Q/I spec- 
trum, the continuum level being rapidly reached as one moves 
away from line center. This explains why quantum interfer- 
ences between HFS F-levels pertaining to different 7-levels, 
whose signatures are negligible in the line core but become 
dominant in the wings, are found to play no significant role 
on the emergent 2/^ profile of this lithum doublet. Moreover, 
we find that such two peaks in the theoretical Q/I profile stand 
up clearly only when the kinetic temperature of the thin atmo- 
spheric region that produces the emergent Q// profile is lower 
than 4000 K. In order to be able to reproduce the width of 
the observed Q/I profile without accounting for non-thermal 
broadening we need at least 6000 K, but for T > 5000 K the 
predicted 2// profile is slightly asymmetric (i.e., with a more 
extended red wing) and only shows the dominant peak due to 
^Li I. 

The outline of the paper is the following. In Section 2, we 
present the atomic model and the structure of the hyperfine 
multiplets, both of ^Li i and of ^Li i. In the same section, we 
briefly present the density matrix formalism which we have 
adopted in order to describe the atomic polarization induced 
in the various levels by anisotropic radiation pumping pro- 
cesses. The optically thin slab model considered in this inves- 
tigation, and the relevant equations describing the polarization 
of the emergent radiation are discussed in Section 3. In Sec- 
tion 4, the resulting theoretical Q/I profile is shown, and the 
physical origin of its spectral features is investigated. Section 
4.1 is dedicated to a discussion of the role of quantum inter- 
ferences between the upper 7-levels of the Di and D2 lines. 
Finally, the sensitivity of the resulting profile to the isotopic 
abundances and to microturbulent magnetic fields is discussed 
in Section 5. Section 6 collects our conclusions with an out- 
look to future research. 

2. THE ATOMIC MODEL 

We adopt a three-level model of Li i consisting in the 
ground level (2s ^Si/2), the upper level of the Di line {2p 
^Pi/2), and the upper level of the D2 line {2p ^P3/2)- The ener- 
gies of the fine structure levels are taken from lRalchenko et~an 
( I2OO8I) . In order to take into account the isotopic efifect, we 
correct the energies of the ^Pi/2 and of the ^P3/2 levels of ^Li i 
using the values of the isotopic shifts in the Di and D2 lines 
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Fig. 2. — Grotrian diagram showing the terms, the fine structure, and the hyperfine structure levels considered in our atomic models for the 'Li and 'Li isotopes 
(splittings are not drawn to scale). On the Grotrian diagrams the HFS components of the Dj and D2 lines of both isotopes are drawn. Panel a shows the laboratory 
positions and relative intensities of all the HFS components, as obtained taking into account the meteoritic relative abundance of the two isotopes. Panels b, c and 
d show in more detail the HFS components of the 'Li i D2 fine {b), the HFS components of the 'Li i Di line together with those of the *Li i D2 line (c), and the 
HFS components of the *Li i Di line id). Note that the scale of the Ijl™'^^ graph is different in the various panels. The zero of the wavelength scale in panels h, c 
and d is taken at 6708 A. 



listed in Table [T] 

The HFS Hamiltonian, describing the interaction between 
the nuclear spin and the electronic angular momentum, can be 
expressed as a series of electric and magnetic multipoles (see, 
for example. iKopfermannI 1 95 8l) . We calculate the energies of 
the HFS F-levels using the values of the magnetic dipole and 
of the electric quadrupole HFS constants (usually indicated 
with the symbols J{ and S, respectively) listed in Table [T] 
We recall that in the absence of magnetic fields, using Dirac's 
notation, the energy eigenvectors can be written in the form 
\aJIFf >, where a represents a set of inner quantum num- 
bers (specifying the configuration and, if the atomic system 
is described by the L-S coupling scheme, the total electronic 
orbital and spin angular momenta), J is the total electronic 
angular momentum quantum number, while F and / are the 
quantum numbers associated with the total angular momen- 
tum operator (electronic plus nuclear: F=J-i-I), and with its 



projection along the quantization axis, respectively. The Gro- 
trian diagrams showing the various HFS F-levels of the two 
isotopes, and the HFS components of the Di and D2 lines are 
shown in the upper panel of Figure |2l In the lower panels of 
Figure |2] the laboratory positions of the various HFS compo- 
nents are shown. Since the isotopic shifts are of the same or- 
der of magnitude as the frequency separation between the two 
D-lines, the Di line of ^Li i falls almost at the same wave- 
length as the D2 line of *Li i (see panels a and c of Figure|2]). 
We observe that in both isotopes the ground level splits into 
two HFS F-levels. The splitting between these two levels is 
much larger than that among the HFS F-levels pertaining to 
the upper levels. As a consequence, the HFS components both 
of the Dl and of the D2 lines can be gathered into two groups. 
Note also that it is not possible to resolve, even in the labo- 
ratory spectra, all the HFS components of the D2 line, since 
their frequency separation is smaller than the natural width 
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TABLE 2 

Wavelengths (in air) and Einstein coefficients of the lines considered; 
mean number of photons and anisotropy factor of the photospheric 
continuum at the wavelengths of the same lines. 



Line 


^(A) 


A(.v-') 




Wy 


Di 


6707.91 


3.72x10' 


0.0115 


0.0992 




6707.76 


3.72x10' 


0.0115 


0.0992 



(compare the magnetic dipole HFS constant of the D2 upper 
level with the Einstein coefficient for spontaneous emission in 
this line, quoted in Table |2l). 

We describe the excitation state of the Li i levels by means 
of the density matrix formalism, a robust theoretical frame- 
work very suitable for treating the atomic polarization (pop- 
ulation unbalances, and quantum interferences among the 
magnetic sublevels) that can be induced, for instance, by an 
anisotropic incident radiation field. In principle, to have a 
complete description of the atomic polarization, one has to 
take into account all the quantum interferences (or coher- 
ences) of the form 



< aJIFf\p\aJ'IF'f' > 



(1) 



where p is the density operator. In this investigation we, first, 
restricted to the 7-diagonal density matrix elements 



< aJIFf\p\aJIF'f' > 



(2) 



or, in other words, we neglected coherences between different 
7-levels. According to LL04, the resulting model atom is re- 
ferred to as multi-level atom with HFS. The statistical equilib- 
rium equations (SEE), and the radiative transfer coefficients 
for a multi-level atom with HFS can be found in LL04. More 
general calculations, that we carried out including the interfer- 
ences between the upper HFS levels of the Di and D2 lines, 
showed that the previous approximation is fully justified for 
the investigation of this lithium doublet. Indeed, as it will 
be shown in Section 4. 1, the theoretical Q// profiles obtained 
through the solution of these more general equations cannot 
be distinguished from those shown in this paper, which, as 
mentioned above, neglect quantum interferences between dif- 
ferent 7-levels. This is due to the weakness of the lithium 
doublet which does not allow for the formation of extended 
wings where such quantum interferences would produce their 
main signatures. 

3. THE OPTICALLY THIN SLAB MODEL 

In order to emphasize the atomic aspects involved in the 
problem, avoiding complications due to radiative transfer ef- 
fects, we consider a horizontal, optically thin slab of Li i ions. 
We assume the slab to be illuminated from below by the ob- 
served photospheric continuum radiation field, that we sup- 
pose to be unpolarized and cylindrically symmetric about the 
local vertical. Under these assumptions, taking a reference 
system with the z-axis (the quantization axis) directed along 
the vertical, only two components of the radiation field tensor, 
through which we describe the incident continuum radiation, 
are non-vanishing: 

n fdQ , fdQ 1 , 

^S(v)=^ — /(v,A() and jI(v)=^ ——{3f/-mv,fi) . 

(3) 

where jj. is the cosine of the heliocentric angle. The former 
quantity is the mean intensity of the incident radiation field 



(averaged over all directions), the second one quantifies its 
degree of anisotropy (unbalance between vertical and hori- 
zontal illumination). Mean intensity and anisotropy degree of 
the radiation field can also be expressed in terms of two non- 
dimensional quantities: the average number of photons per 
mode, n, and the so-called anisotropy factor, w (which varies 
between w = -1/2, for the case of horizontal illumination by 
an unpolarized radiation field without azimuthal dependence, 
and w = 1 for vertical unpolarized illumination). Such quan- 
tities are related to and 7^ by the equations 



n(v) 



2hv^ 



Joiy) , w(v) = V2 



(4) 



We calculate the values of n(v) and w(v) of the photospheric 
continuum at the frequencies of the Li i D-lines, following 
Section 12.3 of LL04 taking h = 0, and using the values of the 
disk-center intensities and of the limb-darkening coefficients 
given by Pierce (2000). The values obtained are listed in Ta- 
ble |2] Once the SEE have been written down and solved nu- 
merically, we can calculate the radiative transfer coefficients. 
We consider the radiation scattered by the slab at 90°, and we 
take the reference direction for positive Q parallel to the slab. 
For the case of a tangential observation, in a weakly polarizing 
atmosphere (e/ » eg, e^/, ey; ?// » riQ,riu,riv,pQ,pu,Pv), 
the polarization of t he emergent radiation is given by (see 
iTruiillo Buenoll2003h 



X(v, il) ^ sx(v, il) i]x(v, n) 
I{v, n) ~ e/(v, n) ~ 77/(v, fi) 



with X^Q,U,V . (5) 



The first term in the right hand side of Eq. ^ represents the 
contribution to the emergent radiation due to selective emis- 
sion processes, the second one is caused by dichroism (selec- 
tive absorption of polarization components). 

The way atomic polarization is distributed among the var- 
ious levels by radiative processes is very similar to the 
case of the sodiurn and barium D-lines, investigatei i by 
Truiillo Bu eno et alj (12002 ) and bv lBeUuzzi et alj ("2007), re- 
spectively. We recall in particular that only the upper level 
of the D2 line can be directly polarized by the anisotropic in- 
cident radiation field. The ground level becomes polarized 
because of a transfer of atomic polarization via spontaneous 
emission in the D2 line, while the Di upper level is polar- 
ized via radiative absorption in the Di line ("repopulation 
pumping"). Indeed, solving the SEE in the absence of mag- 
netic fields, one finds that the D2 line upper level is consid- 
erably more polarized than both the Di-line upper level and 
the ground level. Dichroism therefore can be safely neglected 
as far as D2 is concerned, while it is expected to be more im- 
portant in the Di line. However, given that the polarization 
signal produced by the D 1 line of Li i is expected to be several 
orders of magnitude smaller than the one produced by the D2 
line, we neglect dichroism also in the Di line. We calculate 
therefore the polarization of the emergent radiation through 
the simpler equation 



X{v, £1) sx(v, £1) 
/(v, £1) ~ e/(v, £1) 



(6) 



The emission coefficients that appear in the previous equa- 
tion include only line processes: in order to reproduce the ob- 
served Q/I profile, we need to add the contribution of the con- 
tinuum. Assuming such contributions to be constant across 
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Fig. 3. — Theoretical Q// profile obtained assuming the following values of 
the free parameters: = 200xe™'',eg = The profile is obtained 

for a Doppler width of 60 mA, which corrresponds to T « 3000 K. 

the line we have 

x(v,n) ^ £^(v,n) + 4 

I(v,£l) ~ e/(v, + ' 

where the superscripts "c" and "f ' recall that the correspond- 
ing quantities refer to continuum and line processes, respec- 
tively. The quantities and Sj will be considered as free 
parameters to be adjusted in order to reproduce the observed 
polarization profile. The effect of collisions will be neglected 
throughout this investigation. It should be observed that 
given the extreme weakness of the absorption features of this 
lithium doublet in the solar intensity spectrum, the optically 
thin slab model is expected to be a rather good approximation 
for the modelling of these lines. 

4. THEORETICAL LINEAR POLARIZATION PROFILE OF THE 
EMERGENT RADIATION 

Applying Eq. (Q, and assuming the following values of 
the free parameters Sj - 200 x e"^", where e™^" is the 
maximum value of in the wavelength range considered, 
Sq = 10 X Sj, and a Doppler width of 60 mA, we obtain the 
two-peak linear polarization profile shown in Figure[3] As can 
be observed in the left panel of Figure |4l modifying the con- 
tinuum contribution to the intensity we modify the amplitude 
of the two peaks, without changing the shape of the profile. 

The value that we assumed for the continuum contribution 
(Sj - 200 X e™™) appears to be a good choice, given the 
weakness of the line, and that it allows the peak falling at 
shorter wavelengths (hereafter referred to as the "blue peak") 
to reach the same amplitude (0.02%) as the central peak of the 
Q/I profile observed bv ,Stenflo et al.. (2000.) (see Figure [TJ. 
The value of Sq has been adjusted in order to obtain in the far 
wings the same continuum polarization level as in the obser- 
vation (0.01%). 

Particularly interesting is the sensitivity of the theoretical 
Q/I profile to the value of the Doppler width. The profile 
shown in Figure [3] has been obtained assuming a Doppler 
width of 60 mA, which corresponds, neglecting microturbu- 
lent velocities, to a temperature of 300() K. For the sake of 



simplicity, we consider the same Doppler width for the two 
isotopes, despite their mass difference I4%f. Profiles 
obtained assuming different values of the Doppler width are 
plotted in the right panel of Figure H] Interestingly, we ob- 
serve that the two-peak structure gradually disappears as the 
Doppler width is increased: the observation of a two-peak 
structure could thus provide precise information concerning 
the thermal properties of the thin atmospheric region where 
this weak lithium doublet is formed. 

In the left panel of Figure |5] the same Q/I profile as in 
Figure [3] is plotted together with the profiles expected in the 
hypothetical cases where only ''Li (short-dashed line) or only 
^Li (long-dashed line) were present. The polarization peak 
that is obtained when a single isotope is considered (either 
^Li or ^Li) is, in both cases, due to the corresponding D2 line. 
Indeed, the polarization signals produced by the Di lines, both 
of ^Li I and of ^Li i, are several orders of magnitude smaller, 
and cannot be appreciated on this plot. It is then clear that 
the physical origin of the two-peak structure of the Q/I pro- 
file that we have obtained within our modelling assumption 
lies in the isotopic shift between the two lithium isotopes: the 
two peaks are nothing else but the signals produced by the D2 
lines of ^Li i (blue peak) and of ^Li i (red peak), which fall at 
different wavelengths because of the isotopic shift, and which 
are weighted by the isotopic abundances. 

It is also interesting to observe that the signals due to the 
D2 lines of the two isotopes, as found in the hypothetical 
cases where only ^Li or only ^Li were present, do not have the 
same amplitude (see the left panel of Figure |5j. This is due 
to the fact the two isotopes have different HFS (in particular 
different nuclear spin quantum numbers, and different HFS 
constants, producing different splittings among the various F- 
levels, so that ''Li is less depolarized by HFS than ^Li). As a 
consequence, the relative amplitude of the two peaks that we 
have found in the resulting Q/I profile does not reflect only 
the different abundances of the two lithium isotopes, but also 
their different HFS. Therefore, the presence of HFS in the two 
isotopes has to be taken into account in order to obtain the cor- 
rect relative amplitude between the two peaks. This is clearly 
shown in the right panel of Figure |5] where the profiles ob- 
tained including (solid line) and neglecting (dashed line) HFS 
are plotted. 

4.1. The role of coherences between the upper J -levels of the 
Li I D\ and D2 lines 

Although the Di and D2 lines (both of ^Li i and ^Li i) are 
very close to each other, so that quantum interferences be- 
tween the upper levels of the corresponding transitions are 
not negligible, it is important to point out that, because of 
the weakness of these lines, the continuum level is rapidly 
reached both in the intensity spectrum, and in the Q/I spec- 
trum, and the spectral signatures of such coherences are com- 
pletely lost. This can be easily appreciated in Figure |6l 
where the theoretical 2/^ profiles, obtained both without (up- 
per panel) and with (lower panel) continuum, are plotted ac- 
cording to three different approximations: a) taking into ac- 
count HFS, but neglecting interferences between pairs of HFS 
magnetic sublevels pertaining to different 7-levels (solid line), 
b) neglecting HFS, but taking into account interferences be- 

^ Note that calculations made taking into account the mass difference be- 
tween the two isotopes have not shown any appreciable modification of the 
results. 
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Fig. 4. — Left panel: theoretical Qjl profiles obtained for different values of the continuum {si'). The solid-line profile is the same as in Figure|3] Right panel: 
theoretical Qjl profiles obtained for difterent values of the Doppler width (A/i/j). The solid-line profile is the same as in Figure[3] 
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Fig. 5. — Left panel: theoretical QII profiles obtained {a) considering both isotopes weighted by their meteoritic abundances (solid line), (h) assuming that only 
*Li (100% in abundance) is present (short-dashed line), (c) assuming that only 'Li (100% in abundance) is present (long-dashed line). The vertical dotted lines 
show the wavelength positions of the HFS components of the ^Li i D2 line (first line going from shorter to longer wavelengths), of the ^Li i Di line together with 
the *Li I D2 line (middle line), and of the *Li i D2 line (last line). The solid-line profile is the same as in Figure[3] Right panel: theoretical Qjl profiles obtained 
considering both isotopes, weighted by their meteoritic abundances, taking into account (solid line) and neglecting (dashed line) HFS. The solid-line profile is 
the same as in Figure|3] All profiles have been obtained for a Doppler width of 60 mA. 



tween pairs of magnetic sublevels pertaining to different J- 
levels (dotted line), and c) taking into account HFS, and tak- 
ing also into account interferences between pairs of HFS mag- 
netic sublevels pertaining to different /-levels (dashed line). 
The first approximation is the one that has been applied in 
this paper (see Section 2). A detailed description of the other 
approximations, usually referred to as multi-term atom and 
multi-term atom with HFS c an be found in LL04, and in 
ICasini & Manso Sainzl (l2005h . respectively. 

We start considering the results obtained without contin- 
uum (upper panel). As expected, close to line center there is 
no appreciable difference between the profile obtained includ- 
ing the interferences between the upper levels of the Di and 
D2 lines (dashed-line profile) and the one obtained neglecting 
such interferences (solid-line profile). Actually, we recall that 
such interferences play an important role in the wings of the 
line, being negligible close to the line core. At line center, on 



the other hand, there is an appreciable difference between the 
profiles obtained taking into account (solid- and dashed-line 
profiles) and neglecting (dotted-line profile) HFS. In the far 
wings, on the contrary, as can be observed by comparing the 
dashed- and the dotted-line profiles, there is no difference be- 
tween the results obtained by taking into account or neglecting 
HFS (as known, the effect of HFS vanishes in the far wings 
of the line)^. At these wavelengths, on the other hand, the in- 
terferences between different /-levels become important, and 
their effect can be easily appreciated. In particular, we observe 
that the asymptotic value reached by the profiles obtained tak- 
ing into account such interferences (dashed- and dotted-line 
profiles) is much higher than the one reached by the profile 
calculated neglecting these interferences (solid-line profile). 

^ Note that this is strictly true only if the lower level is not polarized, or 
if it caiTies an amount of atomic polarization much smaller than that of the 
upper level, as in the lines under investigation. 
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Fig. 6. — Upper panel: theoretical fractional polarization profiles obtained 
in the absence of continuum according to three different approximations (see 
the text). The quantities n, w, and AAq have the same values as in Figure[3] 
The vertical dashed lines show the wavelength positions (going from the blue 
to the red) of 'Li i Do, of 'Li i Di and ^Li i D2 (blended), and of *Li i Di . 
Lower panel: same as the upper panel, but including the contribution of the 
continuum. The quantities Ej and Eq have the same values as in Figure [5] 
The profile plotted with solid line is the same as in Figure[3] 

However, if the continuum contributions to Stokes / and Q 
are included, so that a two-peak profile with the same ampli- 
tude as the observed one is obtained, the polarization in the 
wings of the line is rapidly decreased to the continuum level, 
so that the profile obtained by taking into account the inter- 
ferences between different 7-levels cannot be distinguished 
from the one obtained neglecting such interferences (see the 
lower panel of Figure |6] where the solid and dashed lines su- 
perimpose). The only difference that can still be noticed is at 
line center, between the profiles obtained by including and ne- 
glecting HFS (as already shown in the left panel of Figure|5]l. 
We point out that the depolarization of the Q/I profile in the 
wings of the line, which has completely hidden the signatures 
of the interferences between the different 7-levels, is particu- 
larly strong in these lines since a significant continuum con- 
tribution (Sj - 200 X e™^") is needed in order to reproduce 
the amplitude (« 10"^) of the observed Q/I profile. Never- 
theless, as already pointed out, this is a rather realistic value, 
given the extreme weakness of this lithium doublet"*. In con- 
clusion, as anticipated in Section 2, interferences between the 
upper levels of the Di and D2 lines, though present, can be 
safely neglected in the modelling of these lines. 

5. SENSITIVITY TO THE ISOTOPIC ABUNDANCE AND TO A 
MICROTURBULENT MAGNETIC FIELD 

In the left panel of Figure |7] the profiles obtained changing 
the relative abundance of the two isotopes are plotted. We ob- 
serve that as the ^Li abundance is increased, the amplitude of 
the blue peak decreases, while the amplitude of the red one 

"* Note that also for smaller values of £|' the effects of the interferences 
between dift'erent 7-levels turn out to be hidden. 



increases. This behavior suggests the possibility of using the 
relative amplitude of the two peaks to estimate the abundance 
ratio of the two lithium isotopes in the quiet solar atmosphere. 
Note that the two peaks have a slightly different sensitivity to 
the relative isotopic abundance. This is probably due to the 
fact that the blue peak is only caused by the ^Li i D2 line, 
whereas the red peak, though dominated by the ^Li i D2 line, 
is also affected by the ^Li i Di line. Obviously the amplitude 
of the two peaks of this Q/I profile is not sensitive only to 
the isotopic abundance, but also to collisions (here neglected), 
and to the presence of a magnetic field. As it can be observed 
in the right panel of Figure [T] in the presence of a unimodal 
microturbulent and isotropic magnetic field the amplitude of 
the two peaks is reduced because of the Hanle effect. Note 
that the two peaks have a slightly different magnetic sensitiv- 
ity (for example, the relative depolarization that takes place 
going from 1 to 2 G is larger in the red peak). This is due to the 
fact that the two isotopes have different nuclear spin quantum 
numbers and different HFS splittings. The saturation regime 
is reached for magnetic fields of about 50 G. In conclusion, 
a correct estimate of the lithium isotopic abundance from the 
relative amplitude of the two peaks of the Q/I profile would 
require an independent determination of the magnetic field 
strength present in the regions of the solar atmosphere that 
produce the observed scattering polarization. Likewise, for 
a given approximate value of the lithium isotopic abundance 
the two peaks of the Q/I profile of the lithium doublet can 
be used to obtain information on the strength of unresolved, 
hidden magnetic fields in the "quiet" sun. The problem is that 
it is unlikely that such two peaks can indeed be observed (see 
below). 

6. CONCLUDING COMMENTS 

The application of the quantum theory of spectral line po- 
larization leads to the conclusion that only two peaks can be 
expected for the Q/I profile of the Li i resonance doublet at 
6708 A. Interestingly, we find that the strongest peak at about 
6707.75 A is caused by the D2 line of ^Li i, while the weakest 
one at about 6707.9 A is due to the D2 line of ^Li i. Since 
for the lithium doublet in the sun the line opacity is negligible 
with respect to that of the continuum, we find that quantum 
interferences between the HFS levels pertaining to different 
7-levels play no significant role on the emergent Q/I profile. 
Moreover, we have showed that such two peaks in the theo- 
retical 2// profile stand up clearly only when the kinetic tem- 
perature of the thin atmospheric region that produces the main 
contribution to the emergent Q/I profile is sufficiently low 
(e.g., sensibly lower than 4000 K). In order to be able to repro- 
duce the width of the observed Q/I profile, without account- 
ing for non-thermal broadening, we need at least 6000 K, 
but for this kinetic temperature the predicted Q/I profile is 
slightly asymmetric (i.e., with a more extended red wing) and 
shows only the dominant peak due to 'Li i. In this respect, it 
is of interest to note that the height in standard semi-empirical 
models of the solar atmosphere where the continuum optical 
depth at the wavelength of the lithium doublet is unity for a 
line-of-sight with /j = 0.1 is slightly below 200 km (see Fig- 
ure 6 of Trujillo Bueno & Shchukina 2009), where the kinetic 
temperature is about 5000 K. For this reason, we believe that 
the most likely shape of the scattering polarization profile of 
the lithium doublet that the sun can produce should be similar 
to the short-dashed line of the right panel of Figure |4] 

If that is the case (i.e., if the true Q/I profile of the lithium 
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doublet turns out to have only one peak) then it would not 
be possible to apply the "peak-ratio technique" illustrated in 
Figure |2] neither for determining the relative lithium isotopic 
abundances nor for estimating the strength of an unresolved, 
hidden magnetic field in the solar photosphere. Nevertheless, 
the amplitude, FWHM, and asymmetry of the observed Q/I 
profile could still be used for such a purpose through radia- 
tive transfer simulations in a realistic model of the thermal 
and density structure of the quiet solar photosphere. For ex- 
ample, if the relative lithium isotopic abundance is known 
beforehand, the discrepancy between the Q/I profile calcu- 
lated in a three-dimensional hydrodynamical model of the 
solar photosphere and the observed profile could be inter- 
preted in terms of an unresolved magnetic field, assuming 
that the depolarizing impact of elastic collisions is properly 
taken into account in the calculations. Since the relevant 
pumping radiation field here is that of the continuum radia- 
tion, whose anisotropy factor is significantly larger above the 
granule cel l centers than above the intergranular lanes (see 
Figure 2 of Truiillo Buen o et al.ll200 4). the observed Q/I pro- 
file of the lithium doublet (which lacks spatio-temporal reso- 
lution) must be significantly biased towards the granule cell 
centers. Therefore, we anticipate that the inferred magnetic 
field strength would be more representative of the magnetiza- 



tion of the granular plasma. 

Obviously, new and very careful observations of this ex- 
tremely weak polarization signal are urgently needed. In this 
respect, it is of interest to mention that our colleagues from the 
Istituto Ricerche Solari Locarno (IRSOL), Dr M. Bianda and 
Dr R. Ramelli, have akeady initiated such an observational 
program with the Zurich Imaging Polarimeter (ZIMPOL) at- 
tached to the Gregory Coude Telescope of IRSOL. Although 
some of the Q/I profiles they have been able to measure with 
such an excellent facility for high-sensitivity spectropolarime- 
try seem to show a one-peak profile, many more observations 
are needed in order to fully clarify what type of Q/I profiles 
the sun is actually producing in the lithium doublet and in or- 
der to use them for diagnostic purposes. 
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